Full-field quantum correlations of spatially entangled photons 
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Spatially entangled twin photons allow the study of high-dimensional entanglement, and the 
Laguerre-Gauss modes are the most commonly used basis to discretize the single photon mode 
spaces. In this basis, to date only the azimuthal degree of freedom has been investigated ex- 
perimentally due to its fundamental and experimental simplicity. We show that the full spatial 
entanglement is indeed accessible experimentally, i.e., we have found practicable radial detection 
modes with negligible cross correlations. This allows us to demonstrate hybrid azimuthal - radial 
quantum correlations in a Hilbert space with more than 100 dimensions per photon. 



High-dimensional entangled photons are of great interest 
in various areas in quantum information, such as they 
promise high-density encoding of quantum information 
[TJ |2] , are more robust against noise and eavesdroppers 
due to stronger non-classical correlations [3], and in gen- 
eral, present a unique model system for the study of 
high-dimensional entanglement in nature. Entanglement 
in the photon's spatial degrees of freedom is a candi- 
date for this, and it can readily be obtained in the lab 
by spontaneous parametric downconversion (SPDC) of 
an intense laser beam. To explore this high-dimensional 
Hilbert space, we need to discretize this initially contin- 
uous space; due to the paraxial nature of experiments, 
this is usually done using a complete and orthogonal ba- 
sis of transverse optical modes. In quantum information 
with entangled particles, it is crucial that the bipartite 
state shows perfect correlation (or anticorrelation) in the 
used quantum numbers. Paraxial optical modes are also 
required for the implementation of quantum cryptogra- 
phy: they are propagation-invariant and superpositions 
thereof are basically stable. Traditionally, a Gaussian 
basis is employed, in particular the Laguerre- Gaussian 
modes LG p with mode indices £ and p proved to be very 
convenient. These modes factor in an azimuthal phase- 
only part u £ az ((j)) = e li ^ and a radial part u^(r), in 
which the azimuthal part gives rise to the photon orbital 
angular momentum (OAM) of £h |4J. This azimuthal 
part, or OAM entanglement, has been subject to a decade 
of numerous very successful experiments (see, e.g., ref. 
[5j |6]), which is well founded by the fact that most ex- 
perimental setups exhibit rotation symmetry around the 
optical axis. 

The amount of entanglement present in the spatial pho- 
ton pairs can be characterized by the average num- 
ber of entangled optical modes, the Schmidt number 
This Schmidt number K = l/T, k X l is ob- 
tained from the eigenvalues (relative weights) of the 
Schmidt decomposition [9 J of the two-photon field = 
Y, k V^k\uk)s\uk)i, where the \u k ) Sii are the Schmidt 
eigenmodes for the signal or idler photon. Although the 
Schmidt modes have to be calculated numerically in the 
general case, the Schmidt number K can be approxi- 




Figure 1: Schematic experimental setup. The spatially en- 
tangled photons are produced in the nonlinear crystal, whose 
surface is imaged with 7.5 x magnification (fai = 10 cm, 
/z,2 = 75 cm) onto the spatial light modulator (display size 
16 x 12 mm 2 , 800 x 600 pixel). This device is programmed 
to perform the phase modulation required to transform the 
detection mode into the fundamental mode. The far-field of 
the SLM surface is imaged (lOx, 0.2 NA objectives) onto 
the single-mode fiber which is connected to a single pho- 
ton counter. Simultaneous detection events from an entan- 
gled photon pair are selected in a coincidence time window 
of 2.3 ns. The inset shows an exemplary SLM phase pat- 
tern {£ s = 1, p s = 1; £i = — 2, pi = 3) and superimposed the 
magnification-corrected waist of the pump beam (purple) , the 
waist of the detection singlemode fiber (1275 mm, yellow), and 
the detection mode waist (red) w = 1000 |im. 



mated as K = |(^cr+^:) 2 , where b~ x (with b 2 = L\ p /Stt) 
is the phase-matching width, and a is the pump beam 
waist [7J. For our experimental parameters (crystal 
length L = 2 mm, pump beam waist w p = 325 [im, pump 
beam wavelength X p = 413 nm), this number is very 
large: K « 350. However, if only the azimuthal degree of 
freedom is employed (i.e., taking p = [10J) this number 
is significantly lower. We can write the two-photon en- 
tangled state as = Ylt^-oo yf^t%P = tys\-t,P = 0)», 
where \£,p) is a photon with OAM £H and radial quantum 
number p, and the (azimuthal) Schmidt number becomes 
Kclz = 1/ For large K, this can be approximated 



as K az « 2\[K [TTJ[T2]. Direct experimental determina- 
tion of this number has been shown only recently |13| . 
For our case, this number is K az « 37, which is ob- 
viously much lower than the total number of entangled 
modes. The "missing" entanglement becomes accessible 
if also the radial modes are taken into account. There, 
we find a radial Schmidt number (for only one azimuthal 
mode, e.g., for t = 0) K ra d « y/K, which in our example 
is K ra d « 18. The vast majority of the entangled modes 
are radial-azimuthal cross-correlated modes |14j . 
Subject to experimental feasibility, the radial part of the 
Laguerre- Gauss entangled modes is an entanglement re- 
source on equal footing; however, only recently it has 
been investigated in detail theoretically |15| . The LG 
mode functions factor as LG^(r, <fi) = C • u £ az ((j)) • u p r ^ d {r)\ 
the azimuthal part is fully orthogonal in £ and indepen- 
dent on the experimental choice of the detection mode 
waist, and the entangled photons are perfectly anti- 
correlated in I: OAM is conserved in SPDC (In Fig. SI of 
the supplementary information we show that this state- 
ment holds also for higher-order radial modes). There- 
fore, the azimuthal modes are automatically Schmidt 
modes. It turns out that in contrast to this, the radial 
modes do not necessarily represent Schmidt modes |7J 
and we expect to find non-zero quantum correlations of 
detected modes with different p. However, for a proper 
combination of pump-beam, detection-mode and phase 
matching waist w* i = ^4b/a [7, llj, also the radial LG 
modes are Schmidt modes and the cross-correlations for 
p s 7^ pi disappear. In our case, we obtain w* i = 37 fim. 
If we neglect phase matching [15J, w^^O. 
Also experimentally, investigation of the radial correla- 
tions turns out to be more challenging: The spatial cor- 
relations are traditionally investigated using a mode con- 
verter (spatial light modulator or spiral phase plate in the 
case of azimuthal correlations), to transform a certain op- 
tical mode into the fundamental Gaussian, which in turn 
can be tested by sending the photon into a single-mode 
fiber. This works very well for the azimuthal modes, but 
for radial modes, complications occur: for instance, the 
finite acceptance angle of the SMF becomes problematic, 
and there are no perfect spatial modulators which allow 
control over amplitude and phase simultaneously (the or- 
thogonality of p-modes requires amplitude-sensitive de- 
tection) , and careful choices for the detection mode waist 
and the fiber collimator have to be taken [7]. We show 
here that despite of these complications, and under the 
right conditions, the relation between the analyzer modes 
and the pump field becomes nicely visible, and a properly 
correlated mode basis can be obtained. 

Experiment. We generate the spatially entangled pho- 
ton pairs by collinear SPDC in a PPKTP crystal (length 
L = 2 mm) of a LGq laser beam (Kr + , A = 413 nm, beam 
waist at crystal w p = 325 ^m, 50 mW power). 
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Experiment Theory 




5 10 15 20 5 10 15 20 



Figure 2: Quantum correlations between radial modes with 
different p (for t s = ii = 0): Shown are the normalized (di- 
vided by maximum) coincidence count rates (color coded) as a 
function of the radial mode numbers p s (horizontal axis) and 
Pi (vertical axis) of the detection modes. Different rows depict 
results for different detection mode waists as indicated. Left 
column: experimental data, right column: theoretical predic- 
tion. It is clearly visible that the smaller the detection mode 
waist gets, the smaller the off-diagonal counts will be. This is 
a sign that we approach the Schmidt basis for 7—^7*. The 
detection mode waists corresponds to waist ratios of (from 
top to bottom) 7 = 2.4; 3; 4.9, see [15] . 



As sketched in Fig. 1, we image the crystal surface with 
7.5 x magnification using a telescope onto the SLM sur- 
face. The SLM is used under an incident angle of 10 or 
5 degrees; this allows us to use a single SLM for both 
signal and idler photon. The SLM is corrected for phase 
flatness, and we operate not in direct phase modulation, 
but use a blaze towards 2 mrad to further lessen the influ- 
ence of phase errors. The far field of the SLM surface is 
imaged onto the single mode fiber using 10 x objectives, 
with a detect ion- fiber mode waist at the SLM of 1275 [im. 
The fibers are connected to single-photon counters and 
we post-select entangled photon pairs by coincidence de- 
tection (time window 2.3 ns). Since the crystal surface is 
imaged onto the SLM, it is sufficient to discuss the situa- 
tion there. The inset of Fig. 1 shows the resulting waists 
of the pump beam, the detection mode, and the detec- 
tion single-mode fiber; with exemplary phase patterns 
for two different settings of the detection mode quantum 
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numbers. The choice of waists depends on (i) the desired 
ratio 7 = w p /w s ^ which determines the orthogonality 
and overlap with the Schmidt modes (where the ideal 
ratio is 7* = 8.8), (ii) the maximum mode order which 
should be detected, this is connected to the singlemode 
detection fiber mode waist, and (iii) the number of entan- 
gled modes required. Our choice of waists is optimized 
for radial and azimuthal mode numbers up to about 10. 
Our SLM-based mode detectors can not project upon 
perfect LG modes, because the amplitude cannot be 
modulated (this is not possible with conventional SLMs 
|2Q|). This can lead to p- non-orthogonal detection fields 
u l v because / °° arg[LG^ (r) LG l p2 (r)} ^ S PliP2 , and one 
would anticipate that cross-correlations will always ap- 
pear; our results below show that this is not always 
true, and that careful adjustment of the detection mode 
waist allows detection of radially entangled modes with 
negligible cross-correlations. Basically, optical diffrac- 
tion couples phase and amplitude, which helps to obtain 
amplitude-sensitive detection. 

For theoretical calculation of the expected coincidence 
count rate, we apply Klysko's picture of advanced waves 
|16j . The detection field (in the near field of the SLM) 
is determined by the Gaussian amplitude of the single- 
mode detection fiber, and the phase as defined by the 
SLM: Ktl = exp[iarg[LG£] -r 2 /w 2 SMF \. We then 
decompose this field in terms of LG modes: u p f al = 
J2 P > °y ,£ ' LG^ . This expansion contains very high-order 
p-components due to the phase jumps at the zeros of the 
LG polynomial with finite intensity. These singularities 
automatically disappear while weighting the modes with 
their relative weight as produced in SPDC C p ^~ l (Eq. 20 

in [EJ). This results in a' pf l = ^ir^C^J J dr r LG £ p , (r) • 

u xtaii wn i cn allows us to calculate the effective detection 
field at the crystal u p ^ t = J2 p , oy,£ • LG p ,. To obtain 
the coincidence amplitude, we can then simply calculate 
the overlap of the two individual detection fields in the 
signal and idler path: C = f d 2 ru p d s e f s • u p ^ 1 . Finally, 
the experimentally observable coincidence count rate is 
r = \C'\ 2 . We use here the the thin-crystal limit (L — » 0), 
which implies that phase-matching effects are neglected 
(experimentally, we are close to perfect phase matching) . 



Radial-mode correlations. Fig. 2 shows the quantum cor- 
relations of purely radial modes (i s = ii = 0) of down- 
converted photons. We clearly observe, as we decrease 
the detection mode waist, that the off-diagonal elements 
in the correlation matrix decrease. This is expected, 
as mentioned above, the radial cross-correlations disap- 
pear for w Sy i — » w* i . To within our experimental ac- 
curacy, we also reproduce the theoretical results of Mi- 
atto et al. [15] very well. Even minute details of the 
experimental data are reproduced qualitatively well by 
our model (Fig. 2); this suggests that our modelling ap- 
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Figure 3: The influence of detection mode waist: The quantity 
on the y-axis is a measure of how sharply the radial quantum 
correlations peak around p s — p%. For ideal LG mode detec- 
tors, we have W —> 1 for the waist ratio 7—^7*. For small 
waists (w < 500 /im), pixilation effects are non-negligible. For 
larger waists, the agreement between experiment (circles) and 
theory (grey dots) is very good. The top axis indicates the 
ratio 7 of the pump beam waist w p to the detection mode 
waist w = w s ,i of the signal and idler photon. 



proach of expanding the detection field in terms of the 
LG modes provided by the SPDC light, is a sound choice. 
For the case of 500 [im mode waist, we estimate the 
(radial) Schmidt number to be 10.4 (experiment) and 
11.2 (theory). This is less than the expectation men- 
tioned above {K ra d — 18), however, SLM pixilation be- 
comes relevant at such small mode waists. To investigate 
this, we determine a measure of the cross correlations, 
or the width of the diagonal (Fig. 2) around p s = pf. 
w = Y, P r (Ps = Pi = P)/Y, Ps , Pi r (Ps,Pi)- For perfectly 
orthogonal modes, W should be unity. Fig. 3 shows W as 
a function of the beam waist, comparing our theoretical 
simulation with experimental data, again we find good 
agreement. This dependency of the cross-correlations 
on beam waist ratio persists also for higher azimuthal 
modes (£ s = —ii), as shown in Fig. S2 of the supple- 
mentary information. We observe that for a detection 
mode waist smaller than 500 ^m, the "orthogonality" W 
decreases again, this (and the fact that the theoretical 
curve does not reach unity) is a consequence of SLM pix- 
ilation: 500 ^m waist corresponds to ~ 25 pixel of the 
SLM. Our results also demonstrate that the apprehension 
of Miatto et al. [15], that the experimentally accessible 
mode waist ratios 7 are too small, therefore leading to 
strong cross-correlation in p-space, which would imply 
that radial modes are not useful for quantum informa- 
tion, is overcautious: We can adjust the mode waists 
that cross-correlations become negligible. 
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Experiment Theory 
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Figure 4: Quantum correlations between radial modes with 
different p, for given £ (rows) at a fixed waist ratio 7 = 2.4. 
The left column shows experimental data, the right column 
the theoretical prediction. The number of single-particle 
modes involved in these experiments is approximately 100. 



step forward beyond the conventionally used azimuthal 
degree of freedom, or orbital angular momentum entan- 
glement. We find that the radial degree of freedom is 
indeed a useful entanglement resource, if care is taken: 
Our experiments and the theoretical model show that 
the choice of detection mode waist is crucial and has to 
be taken into account carefully; we are able to demon- 
strate the transition to a detection basis where cross- 
correlations disappear, effectively a transition to a quasi- 
Schmidt basis. An important next step will be confirma- 
tion and quantification of the "hybrid" azimuthal-radial 
mode entanglement, which is beyond the scope of this 
paper. If radial modes in spatial entanglement are acces- 
sible, the number of useful entangled modes is roughly 
squared compared to the OAM case, this quadratic in- 
crease in the usable Schmidt number could stimulate new 
experiments like detection- loophole free [T7j Bell tests. 
The higher entanglement density per mode area will also 
enable higher channel capacities in systems where the 
spatial extent is relevant: For the transport of spatially 
entangled photons through optical fibers [18], and also 
through turbulent atmosphere [19J. 
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Radial- azimuthal correlations. Finally, we address the 
question if we can make use of the full azimuthal-radial 
Hilbert space experimentally. We find, independent on 
the radial mode index p, the anticorrelation condition 
£ s = or OAM conservation, is preserved for high 
values (up to ~ 20) of p and £ (see Fig. SI). Fig. 4 
shows p-correlations for a fixed detection mode waist of 
1000 ^m, as in Fig. 2, but this time for different choices 
of i s = —ii = £. Compared to i = 0, we observe negli- 
gible increase of the cross correlations for p s ^ pi, which 
is very encouraging, because this suggests that very high 
dimensional Hilbert spaces become accessible. Addition- 
ally, the p s = Pi correlations get more evenly distributed, 
in agreement with theoretical predictions [15], which also 
increases the number of usable modes. Our results in 
Fig. 4 show two-photon correlations in an approximately 
100 x 100 dimensional Hilbert space. 

Conclusions. In conclusion, we have shown first exper- 
iments with high-dimensionally spatially entangled pho- 
tons in the full Laguerre- Gauss-like basis. We analyze 
the entangled photons in the complete transverse basis 
involving azimuthal and radial correlations; this goes a 
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Supplementary information 

p s = Pi = 5 




SI) Azimuthal quantum correlations at different p s = Pi. We see that the azimuthal quantum correlations are well 
preserved (i.e., that cross-correlations are small), also if the detected state involves high radial quantum numbers. 
The axes in the graphs below are the azimuthal quantum numbers £a (x-axis) and £b (y-axis). The color coded pixels 
represent the coincidence counts, normalized to unity. 
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S2) Degree of correlation in p-space, similar to Fig. 3: Also for higher azimuthal quantum numbers £a = — £b we find 
that, in general, by decreasing the detection mode waist, the degree of correlation is increased. Only for very small 
detection mode waists w < 600 /im we observe again (comp are Fig. 3) an increase of undesired cross-correlations. 



